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We report on a study of diffractive dijet production in pp collisions at <Js = 1.96 TeV using the 
CDF II detector at the Fermilab Tevatron pp collider. A data sample from 310 pb _1 of integrated 
luminosity collected by triggering on a high transverse energy jet, E^ 1 , in coincidence with a recoil 
antiproton detected in a Roman pot spectrometer is used to measure the ratio of single-diffractive 
to inclusive-dijet event rates as a function of x p of the interacting parton in the antiproton, the 
Bjorken-z, x p Bj , and a Q 2 w (E 3 T et ) 2 in the ranges 10 -3 < x p Bj < W 1 and 10 2 < Q 2 < 10 4 GeV 2 , 
respectively. Results are presented for the region of p-momentum-loss fraction 0.03 < £p < 0.09 and 
a four-momentum transfer squared tp > —4 GeV 2 . The tp dependence is measured as a function of 
Q 2 and x p B j and compared with that of inclusive single diffraction dissociation. We find weak x%j 
and Q 2 dependencies in the ratio of single diffractive to inclusive event rates, and no significant Q 2 
dependence in the diffractive tp distributions. 
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We report on measurements performed on diffrac- 
tive dijet production in proton-antiproton collisions at 
y/s = 1.96 TeV at the Fermilab Tevatron in Run II 
using the CDF II detector. The CDF collaboration 
has previously reported several results on low-transverse- 
momentum (soft) and high-transverse-momentum (hard) 
diffractive processes, obtained at the Fermilab Tevatron 
pp collider in Run I (1992-96) [l|-[I3] and in Run II 
(2001-06) [l7l-[2ii|. Three of the Run I results [E H3, [H 
are based on dijet events, J J, produced in single diffrac- 
tion (or single dissociation, SD). These events are char- 
acterized by the presence of two jets in the final state and 
a leading p that escapes the collision intact. A large ra- 
pidity gap, Gp, defined as a region of pseudorapidity 21 1 
devoid of particles, is present between the surviving p 
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and the diffractive cluster, X p , composed of the particles 
from the dissociation of the proton, including J J: 



p + p — > p 



■ Gp + X p , 



X„ 



X + JJ, 



(1) 



where X represents all other particles in the diffractive 
cluster and will be referred below as the underlying event 
(UE). 

Among the results of the Run I diffractive measure- 
ments, the most striking one is the observation of a break- 
down of QCD factorization, expressed as a suppression 
by a factor of O(10) of the diffractive structure function 
(DSF) measured in dijet production relative to that de- 
rived from fits to parton densities measured in diffractive 
deep inelastic scattering (DDIS) at the DESY e-p col- 
lider HERA (see OH). This result was further explored 
by CDF in Run I in studies of other diffractive processes, 
as briefly discussed below. 

In this paper, we study diffractive dijet production 
in Run II with an upgraded CDF detector, incorpo- 
rating special forward detector-components that expand 
the kinematic reach of the variables that define the pro- 
cess. Our goal is to further characterize the properties 
of diffractive dijet production in an effort to decipher 
the QCD nature of the diffractive exchan ge, tradition- 
ally referred to as Pomeron (IP) exchange 22|-[24j|. Be- 
low, we first present a brief physics-oriented summary of 
the diffractive processes studied and the results obtained 
in Run I, which is intended to define the terms and con- 
cepts used in the present study. Then, we restate the aim 
of this study in light of the Run I results, and conclude 
with an outline of the organization of the paper. 

In Run I, several soft and hard diffraction processes 
were studied at y/s = 1800 GeV, and in some cases at 
\fs = 630 GeV. Two types of hard diffraction results 
were obtained: 

(a) gap results: diffractive to nondiffractive (ND) 
cross section ratios using the (large) rapidity 
gap signature to select diffractive events; 

(b) RPS results: diffractive to ND structure func- 
tion ratios using a Roman Pot Spectrometer 
(RPS) to trigger on and measure a leading 
antiproton. 

The diffractive dijet production processes studied are 
shown schematically in Fig. [1] They include (a) sin- 
gle diffraction (or single dissociation, SD), (b) double 
diffraction (or double dissociation, DD), and (c) double 
Pomeron exchange (DPE) or central dissociation (CD). 
For SD results have also been obtained for W |5| , b- 
quark and J/ip [l3[ production. 

Soft diffraction processes studied include SD [2j, 
DD QJ], DPEfor CD) [Hj], and single plus double diffrac- 
tion (SDD) These processes are defined as follows: 

SD pp->p + Gf, + X n 

DD pp-> X p 

DPE or CD pp -)• p 

SDD pp^p + G s 



Gp 
G c 



X n 



X c + G p - 
X c + G v 



Oap 


Jet+Jet 


4> 




Gap 


Jet :: 


4> 


Gap 


Jet+Jet 


Gap 




U 
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FIG. 1: Leading order schematic diagrams and event topolo- 
gies in pseudorapidity (77) vs azimuthal angle (<f>) for diffrac- 
tive dijet production processes studied by CDF: (a) single 
diffraction, (b) double diffraction, and (c) double Pomeron 
exchange; the dot-filled rectangles represent regions where 
particle production occurs. 



where Xp, X p and X c are antiproton-, proton- and 
central-dissociation particle clusters, G c is a central gap, 
and Gx p is a gap adjacent to a proton-dissociation clus- 
ter. 

The results obtained exhibit regularities in normaliza- 
tion and factorization properties that point to the QCD 
character of diffraction [25|. For example, at y/s =1800 
GeV the SD/ND ratios, also referred to as diffractive- 
or gap- fractions, for dijet, W, 6-quark, and J/ip produc- 
tion, as well as the ratio of DD/ND dijet-production, are 
all s» 1%. This value is suppressed relative to standard 
QCD-inspired theoretical expectations (e.g., two-gluon 
exchange) by a factor of O(10), which is comparable to 
the suppression factor observed in soft diffraction rela- 
tive to Regge-theory-based predictions. However, except 
for an overall suppression in normalization, factorization 
approximately holds among various processes at fixed 

In this paper, we report on the dependence on Bjorken- 
x, Q 2 and t, the square of the four-momentum transfer 
between the incoming and outgoing antiproton, of the 
diffractive structure function measured from dijet pro- 
duction in association with a leading antiproton in pp 
collisions at y^i=1.96 TeV and compare our results with 
those from HERA and with theoretical expectations. 

The paper is organized as follows: in Sec. [H] we discuss 
the method we use to extract the diffractive structure 
function; in Sec. IIIII we describe the experimental appa- 
ratus; in Sees. [TV] and [V] we discuss the data sets and 
data analysis; in Sec. [VTJ we present the results; and in 
Sec. IVIIIwe summarize the results and draw conclusions. 



II. METHOD 



The cross section for inclusive dijet production in pp 
collisions can be written as 



X r , 



dxp~ dx rp dt 



F™\x p ,Q 2 ) F™\x p ,Q 2 ) 



do 



j j 



.In 



dt 



(2) 
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where Xp (x p ) is the Bjorken-x, defined as the forward- 
momentum fraction of the interacting parton of the an- 
tiproton (proton), Q 2 is the factorization hard scale, 
F^ cl (x p , Q 2 ) and FjJ cl (x p , Q 2 ) are structure functions, 
and djj is the scattering cross section of the two partons 
producing the final-state jets. The Ujj depends on Xp 
and Xp, and also on t, the square of the four-momentum 
transfer between the interacting partons, on Q 2 , and 
on a s (Q 2 ), the strong-interaction running coupling con- 
stant. 

The structure function relevant for dijet production is 
a color-weighted combination of gluon (g) and quark (q) 
terms given by: 



F^f(x,Q 2 ) = x 



g{x,Q 2 



9 ^ 



<li(x,Q 2 ) 



(3) 



where x refers to In analogy with Eq. @, the 

differential cross section for diffractive dijet production 
can be written as 25] 



F?»(x p ,Q 2 ,t,t) F;f(x p ,Q 2 ) da. 



dxpdxpdid^dt 



dt 



(4) 

where Fjj is the diffractive structure function, which 
in addition to the usual dependence on Xp and Q 2 also 
depends on £, the fractional forward momentum loss of 
the antiproton, and on t. 

The jet energies measured in the CDF detector must 
be corrected for various detector effects, which depend 
on the jet energy and rj-<p coordinates due to differences 
in calorimeter subsystem designs and calorimeter inter- 
faces [13 ■ The correction factor generally increases as 
the jet energy decreases. To avoid systematic uncertain- 
ties associated with estimating corrections using Monte 
Carlo simulations, particularly for diffractively produced 
jets of relatively low Et, we measure ratios of SD to in- 
clusive production as a function of Bjorken-x and Q 2 , 
for which jet-energy corrections due to detector effects 
nearly cancel out. 

Since the single-diffractive cross section is a small frac- 
tion of the inclusive one, iS 1% of cr mcl , we refer below to 
"inclusive" and "nondiffractive" (ND) dijet-production 
interchangeably. Furthermore, we also use this approx- 
imation in classifying event samples and associated pa- 
rameters, as for example between F^ cl and Fjj - 

The diffractive structure function is obtained by mul- 
tiplying the ratio Rsd/nd^i^) 01 the SD to ND event 
densities, n^f(x,Q 2 ^,t) and nfP(x,Q 2 ), by the ND 
structure function, Fj^ D (x, Q 2 ): 



F^(x,Q 2 ,C,t) = R S d/nd(z,£,«) x F" D (x,Q 2 ). (5) 

This method of measuring diffractive structure functions 
relies on the LO QCD expectation that cross sections are 
proportional to structure functions: 



Rsd/ndO^Q 2 :^) 



(x,Q 2 ,^t) FfP(x,Q 2 ,^t) 



n» D (x,Q 2 ) F^(x,Q 2 ) 



j j 



(6) 

Next-to- leading-order corrections to F^ D (x, Q 2 , £, t) 
obtained by this method are expected to be of 

0(10%) H3. 



III. EXPERIMENTAL APPARATUS 

The CDF II detector is equipped with special forward 
detectors [H, [28l - [3l| designed to enhance the capabili- 
ties for studies of diffractive physics. These detectors 
include the RPS, the Beam Shower Counters (BSC), and 
the MiniPlug (MP) calorimeters. The RPS is a scin- 
tillator fiber tracker used to detect leading antiprotons; 
the BSC are scintillator counters installed around the 
beam-pipe at three (four) locations along the p (p) direc- 
tion and are used to identify rapidity gaps in the region 
5.5 < |?7| < 7.5; and the MP calorimeters [3(| are two 
lead-scintillator based forward calorimeters covering the 
pseudorapidity region 3.6 < < 5.1. The forward de- 
tectors include a system of Cherenkov Luminosity Coun- 
ters (CLC) [32j . whose primary function is to measure the 
number of inelastic pp collisions per beam-bunch crossing 
and thereby the luminosity. The CLC covers the range 
3.7 < < 4.7, which substantially overlaps the MP 
coverage. In this analysis, the CLC is used for diagnos- 
tic purposes, and also to refine the rapidity-gap definition 
by detecting charged particles that might penetrate a MP 
without interacting, yielding a pulse-height smaller than 
the MP tower thresholds. 

Figure [2] shows a schematic plan view of the beamlinc 
elements and forward detectors along the outgoing an- 
tiproton beam direction. The RPS comprises three Ro- 
man pot stations located at a distance of ~ 57 m from 
the nominal interaction point (IP) along the outgoing p 
direction. Each station is equipped with one Roman Pot 
Trigger (RPT) counter and an 64-channel scintillator- 
fiber tracker, which can be used to reconstruct tracks 
both in the X (horizontal) and Y (vertical) coordinates. 
The tracking information, in conjunction with the inter- 
action point coordinates, provided by the central detec- 
tor, can be used to calculate the variables £ and t of 
the recoil antiproton. During running, the RPS detec- 
tors were brought up to a distance of X sa 1 cm from 
the outgoing antiproton beam, with the (more intense) 
incoming proton beam positioned ~ 2 mm farther away. 

In Run I, in which data were collected at lower lumi- 
nosities, the RPS detectors were positioned at approx- 
imately the same distance from the antiproton beam, 
but the more intense proton beam was in between the 
antiproton beam and the RPS detectors. Background 
from beam-halo particles from the proton beam was lim- 
iting the RPS acceptance at low t. Reversing the polar- 
ity of the electrostatic beam separators (EBS) in Run II, 
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FIG. 2: Schematic plan view of the CDF II central (CDF-II) and forward detectors (acronyms are defined in the text). 



brought the antiproton beam in between the proton beam 
and the RPS. This enabled data taking at the higher in- 
stantaneous luminosities of Run II at y/s = 1.96 TeV 
with approximately the same t acceptance as in Run I at 
1.8 TeV. 

The CDF-I and CDF-II main detectors are multipur- 
pose detectors described in detail in Refs. [HI, [34j . The 
most relevant components of CDF II for this analysis 
are the charged particle tracking system and the central 
and plug calorimeters. The tracking system consists of 
a silicon vertex detector surrounded by the central outer 
tracker (COT) [35j . The part of the silicon vertex detecor 
used is the SVX II (36|, composed of double-sided mi- 
crostrip silicon sensors arranged in five cylindrical shells 
of radii between 2.5 and 10.6 cm. The COT is an open- 
cell drift chamber consisting of 96 layers organized in 
8 superlayers with alternating structures of axial and 
±2° stereo readout within a radial range between 40 and 
137 cm. Surrounding the tracking detectors is a super- 
conducting solenoid, which provides a magnetic field of 
1.4 T. Calorimeters located outside the solenoid arephys- 
ically divided into a central calorimeter (CCAL) [371 [38|. 
covering the pseudorapidity range \r]\ < 1.1, and a plug 
calorimeter (PCAL) [3{|, covering the region 1.1 < |ry| < 
3.6. These calorimeters are segmented into projective 
towers pointing back to the central interaction point, 
with granularity Ar] x A</> « 0.1 x 15°. 



IV. DATA SAMPLES AND EVENT SELECTION 

The Run I diffractive dijet results were obtained from 
data samples collected at low instantaneous luminosities 
to minimize background from pileup events. Constrained 
by statistics, the Q 2 dependence of the DSF was mea- 
sured over a limited range. Moreover, due to uncertain- 
ties in the beam position at the RPS location, it was 
difficult to reliably extract normalized ^-distributions. In 
Run II, our goal was to obtain high statistics diffractive 
data samples with low pileup backgrounds from which to 
extract the Q 2 and t dependence over a wide range. To 
achieve this goal, we built special forward detectors, de- 



scribed above in Sec. lIIIl implemented dedicated triggers, 
and developed analysis techniques for background sub- 
traction and RPS alignment. In this Section, we present 
the data samples and event selection requirements ap- 
plied to produce the data sets from which the results 
were extracted. 

This analysis is based on data from an integrated lu- 
minosity L fn 310 pb" 1 collected in 2002-2003. Only 
"good runs" with a minimum L sa 10 nb _1 are used, se- 
lected based on beam conditions, detector performance 
criteria, and the requirement that the following detec- 
tor components be functional: CCAL, PCAL, MP, CLC, 
and BSC. On-line event selection is accomplished with 
a three-level trigger system (Ll, L2, L3) which accepts 
soft-interaction events as well as hard interaction events 
containing high Et jets. The latter are selected at the 
trigger level by requiring at least one calorimeter tower 
with E T > 5, 20, or 50 GeV within \rj\ < 3.5. Leading 
antiprotons with fractional momentum-loss in the range 
0-03 < % < 0.09 were triggered on using a three- fold 
coincidence of the RPT counters. 

Jets were reconstructed using a cone algorithm, in 
which the transverse energy of a jet is defined as E T et = 
Y>iEiSm(9i) with the sum carried over all calorimeter 
towers at polar angles 9i within the jet cone. The mid- 
point algorithm [40( was used, which is an improved it- 
erative cone clustering algorithm based on calorimeter 
towers with Et > 100 MeV. The jets were corrected for 
detector effects and for contributions from the underlying 
event (UE) [26(. Dijet candidate events were required to 
have at least two jets with Et > 5, 20, or 50 GeV de- 
pending on the event sample, and \rj\ < 2.5. 

In order to explore the region of large-transverse- 
energy jets of relatively low cross sections, data samples 
of RPS-triggered events in conjunction with the pres- 
ence of a jet with E T et >5, 20, or 50 GeV in CCAL or 
PCAL were also studied. These samples are referred to 
as RPS-Jet5, RPS-Jet20, and RPS-Jet50. Correspond- 
ing ND samples (Jet5, Jet20, and Jet50) were used for 
comparison. 

The majority of the data used in this analysis were 
recorded without RPS tracking information. For these 
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data, the value of £ was evaluated from calorimeter in- 
formation and will be referred to as £,p AL (see Sec. lVCp . 
The £,p' AL was then calibrated against £ obtained from 
the RPS, £,p PS , using data from runs in which RPS track- 
ing was available (see Sec. IV Ep . Below, we list the def- 
initions of the triggers used in data acquisition and the 
data selection requirements applied in obtaining the data 
samples for this analysis. 

The following trigger definitions are used: 

• RPS: triple coincidence among the three RPS trig- 
ger counters in time with a p gate; 

• RPStrack: RPS with RPS tracking available (in- 
cluded in the RPS trigger); 

• J5 (J20, J50): jet with E J T et > 5 (20, 50) GeV in 
CCAL or PCAL; 

• RPS-Jet5 (Jet20, Jet50): RPS in coincidence with 
J5 (J20, J50). 

The data selection requirements are listed below: 

(i) good-run events: accepts events from runs with no 
problems caused by hardware or software failures 
during data acquisition; 

(ii) flr significance [4l|: selects events with missing 
transverse energy significance SV, < 2 to reject jet 
events in which there is pr due to energy loss in 
calorimeter cracks and/or events with jets and (un- 
detected) neutrinos, such as from IT' -¥ Iv + jets; 

(iii) N(jet) > 2: accepts events with two jets of > 5 
GeV within \r]i et \ < 2.5; 

(iv) splash veto: rejects events that cause splashes (large 
number of hits) in the RPT counters; 

(v) RPT: rejects events that are triggered by accidental 
(not due to the traversal of a single particle) RPS 
counter coincidences (less than 0.1%); 

(vi) SD (0.03 < ^ AL < 0.09): accepts SD events with 
good efficiency, while rejecting backgrounds from 
pileup events consisting of a soft SD event that trig- 
gers the RPS and a ND dijet event. 

Table U lists the number of events surviving these require- 
ments when applied successively to the data. 

The splash events were studied using the RPS tra ck 
data sample before applying the splash veto require- 
ment. The sum of the ADC counts of the three RPT 
counters, SumRPT = £f =1 RPT ADC , is used to re- 
ject splash events. Figure [3ja) shows SumRPT for 
events with/without a reconstructed RPS track. In the 
"RPS track" histogram the peak at ~ 3000 ADC counts 
is attributed to a single minimum-ionizing particle (MIP) 
traversing all three RPT counters. 2-MIP and 3-MIP 
shoulders are also discernible at ~ 6000 and ~ 9000 ADC 



counts, respectively. Events labeled as "No RPS track" 
are the splash events. These events dominate the re- 
gion of SumRPT > 5000. Detailed studies using various 
event samples indicate that splash events are likely to be 
due to high-£p diffractive events for which the p does not 
reach the RPS, but rather interacts with the material of 
the beam pipe in the vicinity of the RPS producing a 
spray of particles causing the splash. In the region of 
SumRPT < 5000, the "No RPS track" distribution has 
a peak similar to that seen in the events with a recon- 
structed RPS track. These events are interpreted as good 
events for which the track was not reconstructed due to 
either malfunction or inefficiency of the fiber tracker. 

Events with SumRPT> 5000 are rejected. The re- 
tained events contain approximately 77% of the RPT 
sample, as can be qualitatively seen in Fig. EJb) . These 
events, after applying the RPT selection requirement, 
constitute the SD event samples listed in Table |TJ Any 
possible inefficiency caused by the SumRPT cut is taken 
into account by folding a 6% uncertainty into that of the 
extracted cross section (see Sec. lVIB"T| Table InT]) . 

V. DATA ANALYSIS 

Analysis details are presented in seven subsections or- 
ganized by measurement topic, as follows: 

A - Dynamic alignment of Roman pos spectrom- 
eter 

B - Trigger efficiency of Roman pot spectrometer 

C - Antiproton momentum-loss measurement 

D - Multiple interactions 

E - Calibration of £? AL 

F - MiniPlug contribution to £<? AL 

G - Beam-halo background 

A. Dynamic alignment of Roman pot spectrometer 

The values of both £ and t can be accurately deter- 
mined from RPS reconstructed track coordinates and the 
position of the event vertex at the Interaction Point using 
the beam-transport matrix between the Interaction Point 
and RPS. Crucial for this determination is the detector 
X-Y alignment with respect to the beam. Below, we de- 
scribe a method developed to dynamically determine the 
alignment of the RPS detectors during the RPStrack data 
collection period. As described in Sec. IV El the resulting 
^? PS distribution is used to calibrate AL . 

The dynamic alignment method is illustrated in Fig. [4j 
where the curve represents a fit to the data (after align- 
ment) with a form composed of two exponential terms, 

^ = N norrn (A ie blt + A 2 e b ^) , (7) 
where N norm is an overall normalization factor. 
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TABLE I: The number of events surviving successive selection requirements by data sample. 



Selection requirement 



RPS 



RPS-Jet5 



RPS-Jet20 



RPS-Jet50 



Trigger level 
Good-run events 

fir significance: S> ee E t /- s /T.E^, < 2 
N{jet) > 2: Elf > 5 GeV, |V' 2 j < 2.5 
Splash veto 
RPT 

SD (0.03 < ^ AL < 0.09) 



1 634 723 
1 431 460 
1 431 253 
59 157 
27 686 
27 680 
1 458 



1 124 243 
955 006 
950 776 
557 615 
259 186 
259 169 
20 602 



1 693 644 
1 421 350 
1 410 780 
1 168 881 
541 031 
541 003 
26 559 



757 731 
561 878 
539 957 
521 645 
215 975 
215 974 
4 432 



Alignment is achieved by seeking a maximum of the 
da/dt distribution at t = 0. The implementation of the 
alignment method consist of introducing software offsets 
Xoffset and Y ff S et in the X and Y coordinates of the RPS 
detectors with respect to the physical beamline, and it- 
eratively adjusting them until a maximum for da/dt at 
t = (or of the dominant slope b-y) is found at the (X,Y) 
position where the RPS fiber tracker is correctly aligned. 
Results for such a fit are shown in Fig.[SJ The accuracy in 
AX and AY of the RPS alignment calibration obtained 
using the inclusive data sample, estimated from Gaussian 
fits to the distributions in Fig. [5] around their respective 
minimum values, is ±60 /im. This is limited only by the 
size of the data sample and the variations of the beam 
position during data taking. The contribution of the lat- 
ter is automatically folded within the overall uncertainty 
as determined by this method. 



B. Trigger efficiency of Roman pot spectrometer 

Because of radiation damage, the position of the min- 
imum ionizing particle (MIP) peak of the three RPT 
counters shifts toward smaller ADC values as the inte- 
grated luminosity, L, of the data sample increases. The 
same behavior is observed in all three trigger counters. 
As a direct consequence, the efficiency of the RPT triple 
coincidence decreases as a function of L. The RPT ef- 
ficiency, e rpt! is measured using a sample of minimum- 
bias (MB) data from 175.6 pb _1 of integrated luminosity 
collected with a CLC p -CLCp coincidence trigger. The 
ADC distribution of each RPT counter was determined 
for various periods of data-taking by triggering with the 
other two RPT counters. Results are shown in Table |TT] 
for nine data sets of approximately equal integrated lu- 
minosity obtained by subdividing the MB data sample. 
We evaluate crpt from the number of events with at 
least 1000 ADC counts in each of the three RPT coun- 
ters (the trigger requirement) divided by the number of 
events with at least 500 ADC counts in each counter, the 
lowest ADC value found among MIP peaks in all three 
counters and in all data sets: 



e RPT = JVi. 2 . 3 (ADC > 1000)/JVi. 2 . 3 (ADC > 500). (8) 



TABLE II: Data sets of MB events collected with a 
CLCp-CLCp coincidence requirement at various periods dur- 
ing the data-taking run, corresponding integrated luminosi- 
ties, L, and triple-coincidence RPT counter efficiencies, erpt- 



Data set 


L (pb- 1 ) 


ERPT 


set 


12.9 


0.78 ± 0.08 


set 1 


24.0 


0.75 ±0.08 


set 2 


20.3 


0.69 ±0.07 


set 3 


6.4 


0.57 ±0.06 


set 4 


29.2 


0.51 ±0.05 


set 5 


16.3 


0.46 ±0.05 


set 6 


18.9 


0.48 ±0.05 


set 7 


25.5 


0.43 ±0.04 


set 8 


22.1 


0.40 ±0.04 



An uncertainty of 10% is assigned to €rpt to account 
for variations due to the choice of the lowest ADC value 
of the MIP peak as determined from an analysis of the all 
ADC distributions. The degradation of the RPT coun- 
ters is taken into account by dividing the number of ob- 
served SD events by the RPT efficiency corresponding to 
the data-taking period to correct for RPT inefficiencies. 

C. Antiproton momentum loss measurement 

The momentum loss of the p, £^ AL , is calculated using 
the pseudorapidity rf and transverse energy E % T of all the 
towers of the CCAL, PCAL, and MP calorimeters, 

1 AW 

^ CAL = 7i E ^e-"\ (9) 

* i— 1 

where the sum is carried out over towers with Et > 
100 MeV for CCAL and PCAL and E T > 20 MeV for 
MP. Calibration issues of £^ p AL are addressed in Sec. IV El 
The resulting £^? AL distributions are shown in Fig. [6j The 
filled circles in this figure represent the £^ AL distribution 
of the entire RPS-Jet5 data sample after applying all the 
selection requirements listed in Table [IJ This distribution 
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FIG. 3: SumRPT distributions for events of the SD RPS tra ck 
data sample, (a) Number of events and (b) fraction of events 
as a function of the sum of ADC counts of the three Ro- 
man Pot Trigger (RPT) counters, SumRPT. In (a), the peak 
at SumRPT~ 3000, attributed to a minimum ionizing par- 
ticle (MIP) signal, is prominent in the "RPS track" entries, 
where two-MIP (three-MIP) shoulders at ~ 6000 (9000) ADC 
counts are also discernible. The "No RPS track" entries show 
a smaller MIP peak and a broad plateau, which is attributed 
to "splash events" (see text). 
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FIG. 4: t distributions for SD RPStrack events of recon- 
structed RPS tracks in the range < t < 1 GeV 2 for the 
nominal X position and for positive X g se t shifts relative to 
the beamline. 
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shows two peaks: one in the region of AL < 10 1 la- 
beled SD, and a broader one at £? AL > 10 _1 labeled 
BG (background). The data points represented by filled 
squares are from the inclusive Jet5 data sample normal- 
ized (rescaled) to the number of the RPS-Jet5 BG events 
in the region 3 x 1CP 1 < £,p AL < 2, and those repre- 
sented by filled stars are from RPSt ra ck'Jet5 data sam- 
ple events with a reconstructed RPS track normalized 
(rescaled) to the number of RPS-Jet5 events in the SD 
region of 3 x 1(T 2 < £? AL < 9 x 1CT 2 . 

The RPS-Jet5 data distribution is interpreted as fol- 
lows: 

a. Single diffractive region. The events in the SD 
region are mainly due to SD dijet production. The de- 
clining trend of the distribution below £^ AL ~ 5 x 10~ 2 
occurs in the region where the RPS acceptance is de- 
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FIG. 5: The value of —bi of the t-distribution parametrization 
of Eq. (O normalized to its minimum value vs (a) Y fr sc t and 
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FIG. 6: The calorimeter-measured antiproton-momentum- 
loss fraction, £^ AL , for three data samples: (a) RPS-Jet5 
(filled circles), (b) RPStrack- Jet5-rescaled (filled stars), and 
(c) Jet5-rescaled (filled squares) . The regions labeled SD and 
BG are respectively dominated by SD with no pileup and 
background events consisting of a ND dijet event and a SD 
event that provides the RPS trigger. The RPStrack-Jet5 and 
Jet5 distributions are normalized to the RPS- Jet5 distribution 
in the SD and BG regions, respectively. 



creasing. 

b. Background region. The events in the BG region, 
which form the broad peak centered at £^ AL ~ 0.5, are 
mainly due to overlaps of a ND dijet event and a soft 
diffractive interaction that triggered the RPS but yielded 
no reconstructed vertex. A negligibly small fraction of 
the events in this region are due to SD dijet events su- 
perimposed with at least one MB event. Both types of 
overlap events should yield a value of £^ AL ~ 1, but the 
expected nearly 5-function distribution is smeared by the 
resolution of the energy measurement in the calorime- 
ters and shifted toward lower £? AL -values by particles 
escaping detection either in the areas of the calorimeter 
interfaces or due to the imposed energy thresholds. 

As the rate of overlap events increases with instanta- 
neous luminosity £, the ratio of BG/SD events in the 
£? AL distribution of the RPS-Jet5 event sample is ex- 
pected to decrease with decreasing C. We verified this 
effect in the RPS-Jet5 data when binned in intervals of 
different £. 

c. Diffractive events. The excess of the RPS-Jet5 
events over the Jet5 distribution in the SD region is 
mainly due to diffractive production with no pileup. The 
fraction of ND events in the region of 0.03 < C^ AL < 0.09 
is ss 12%. As the RPS acceptance depends on both £p 
and t, and t is not measured for each event, the back- 
ground from ND events in the SD region is accounted 
for in a simple MC simulation designed to calculate the 
acceptance on an event-by-event basis. Inputs to this 
simulation are the £ and t distributions measured from 



data in which RPS tracking is available. Each event gen- 
erated within a given bin of £ RPS is weighted by a factor 
equal to the ratio of the total number of events over the 
number of SD signal events in the corresponding ^ AL 
bin of the RPS-Jet5 data plotted in Fig. [6] 



D. Multiple interactions 

Effects of multiple interactions in the same bunch 
crossing are handled by using the fraction of one- 
interaction events in the data sample, namely the fraction 
of events originating from bunch crossings that have just 
one pp interaction, to normalize the SD/ND event ra- 
tio in a multiple interaction environment. This fraction 
is estimated from a PYTHiA-Monte-Carlo-generated sam- 
ple of events containing the appropriate run-dependent 
fraction of multiple interactions for each run. The in- 
stantaneous luminosities of the various subsamples com- 
prising this dataset vary within the range ~ 2 < C < 
4 x 10 31 cm _2 s _1 . The ratio of the number of one- 
interaction events, TV 1 - 111 ', to all events, 7V a11 , decreases 
with increasing C following a distribution well described 
by an exponential expression 



N 



exp[-c-£], 



(10) 



where c = (0.34 ± 0.06) x 10" 

An offline analysis vertex requirement (cut) accepting 
events with only one reconstructed vertex would elimi- 
nate a large fraction of ND overlap events, but might also 
reject diffractive events due to vertex reconstruction inef- 
ficiencies, which depend on the event activity and there- 
fore on both £ and the number of overlapping events. To 
avoid biasing the £^ AL distribution, we apply no vertex 
cut and correct for the ND overlap event contamination 
when we evaluate cross sections. 



E. Calibration of Q 



The measurement of £p 



ip AL using Eq. © is calibrated by 
comparing with the RPS measurement of £ PPS . The com- 
parison is made using the RPStrack subsample of events 
for which tracking information is available. The ^? PS 
distribution for these events is shown in Fig. [7] 

Figure [5] shows a two-dimensional scatter plot of AL 
vs £? PS for the events with a reconstructed RPS track. 
An approximately linear relationship between £^ AL and 
^ is observed along the peak of the distribution of en- 
tries in the region of £? AL < 0.1 and 0.03 < <^ PS < 0.09, 
where the diffractive events are expected. As already 
discussed, the entries with £^ AL > 0.1 are mainly due 
to a ND dijet event with a superimposed soft SD over- 
lap event. Without the measurement of £^ AL the over- 
lap events would contribute a large and irreducible back- 
ground to the SD events. 
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FIG. 8: Scatter plot of £§ AL vs ^ PS for all RPStrack data. 



For a quantitative calibration of £^ AL the data are di- 
vided into bins of width A^ PS = 0.005 and the £? AL 
values in each bin are fit with a Gaussian distribution. 
Figure M shows such a fit for the bin of 0.055 < < 
0.060. The ratio of the half-width to the average value 
of the fit curve is found to have the constant value of 
S£,p Ah /£,p AL « 0.3 over the entire £? AL range of this 
data-sample. 

The results of £° AL vs £? PS for all the A^ AL bins 
within the plateau region of 0.045 < ^ PS < 0.09 of the 
£p AL distribution of Fig.[7]are plotted in Fig. QUI A linear 



relationship is observed of the form Q 



CAL _ „o 



RPS 



with p° = 0.007 ± 0.002 and p 1 = 0.97 ± 0.04. The small 
size of the deviations of p° and p 1 from the ideal values of 
p° = and p 1 — 1 is testimony to the success of the dy- 
namic alignment of the RPS, which affects ^ PS , and the 
calorimeter energy calibrations and threshold settings, 



FIG. 10: Median values of £^ AL obtained from fits to RPStrack 
data in different ^ AL bins of width A^ AL = 0.05 (horizontal 
bars) fitted linearly in the region 0.045 < £^ AL < 0.09. 



which influence AL . These deviations do not affect the 
referral of the measurement of £,p AL to that of £|^ PS to 
obtain an RPS-based value of t; p in the data for which no 
RPS tracking is available. 



F. MiniPlug contribution to Q 



CAL 



The measured value of & receives contributions 
from all the CDF calorimeter subsystems. The largest 
contribution at large £p comes from PCALp and MPp, 
the PCAL and MiniPlug calorimeters on the p side. The 
contribution of the MPp is particularly important for re- 
ducing the background due to ND dijet events overlapped 
by a soft SD recoil p detected in the RPS. Figure [TT] shows 
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FIG. 11: The ^ AL distribution for RPS-Jet5 and Jet5 data 
evaluated without the MiniPlug calorimeter contributions. 
The diffractive peak is less discernible in the RPS-Jet5 dis- 
tribution, and the Jet5 background is larger than in Fig. [5] 
where the MiniPlugs are included in the evaluation of £^ AL . 



FIG. 12: The ^^cal distribution in zero-crossing (ZC) events 
(circles), ZC with no reconstructed vertices (squares), and ZC 
with a RPS trigger (triangles). 



the £,p distribution for the RPS-Jet5 and Jet5 data 
without the contribution from the MP calorimeters. In 
comparison with Fig. [SJ the BG peak is shifted towards 
the SD region and contributes a considerably larger ND 
background to the SD events. 

The MP calibration was performed by comparing Et 
distributions between data and Monte Carlo generated 
events. The slopes of the dN/dEr distributions of the 
data were first equalized among the MP towers within 
a given polar angular range (ry-ring), and then normal- 
ized to the corresponding slopes obtained from Monte 
Carlo generated events in the same 77-ring. The reso- 
lution in the measurement of Et was estimated to be 
SEt = ±30% using data obtained with a MP prototype 
calorimeter exposed to high energy positron and pion 
beams [28| . Varying the MP energy by ± 30% yields 
a shift of Alog 10 £? AL = ±0.1, which is comparable to 
the bin width used in the analysis (see the data plotted 
in Fig. [B]). 

G. Beam- halo background 

The background in the £^ AL distribution due to beam- 
halo (particles at relatively large distances off the beam 
axis) was studied using a trigger provided by the Teva- 
tron clock at the nominal time the p and p beam bunches 
cross the z = position at the center of CDF II regard- 
less of whether or not there is a pp interaction. This 
trigger, which is referred to as "zero-crossing," leads to 
various data samples depending on additional conditions 
that may be imposed. 

Three zero-crossing event samples were analyzed: 

(a) zero-crossing inclusive; 



(b) zero-crossing events with no reconstructed vertex; 

(c) zero-crossing events with a RPS trigger. 

Figure shows the £^ AL distribution for these sam- 
ples. Three regions of £^ AL are of interest: 

CO £p AL > 0.1, where ND events dominate; 

(ii) 10 -3 < £^ AL < 0.1, where SD events are expected; 

(hi) £,p AL < 10~ 3 , where "empty" events with one to a 
few CCAL or PCAL towers above threshold due to 
beam- halo particles and/or due to calorimeter noise 
may contribute to £^ AL . 

The contribution of a single CCAL or PCAL tower 
at 77 = with Et — 0.2 GeV (the threshold used) is 
estimated from Eq. [3] to be £p AL = 1 X 10~ 4 , which is 
well below the £-range of the RPS acceptance. At £j? AL ~ 
3 x 10~ 4 , we observe a background peak corresponding 
to an average of ~ 3 CCAL towers at an event rate of 
~ 5% of the diffractive signal concentrated at £^ AL ~ 
5 x 10 ~ 2 . We estimate that an upward fluctuation by 
a factor of 100 would be required for this background 
to compete with the overlap background already present 
within the diffractive region of 3 x 10~ 2 < AL < 9 x 
10~ 2 . Because of the negligible probability of such a 
fluctuation, no correction is applied to the data for beam- 
halo background. 



VI. RESULTS 

Our results are presented in three subsections IVI A[ 
IVI CI In IVI A\ we discuss certain kinematic distributions 
that establish the diffractive nature of the events in our 
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FIG. 13: The distribution of the mean dijet transverse energy FIG. 15: The distribution of the azimuthal angle difference 
Et = {E 1 ^ 1 + £f* )/2 for SD and ND events. The ND between the two highest E T jets for SD and ND events, 

distribution is normalized to the total number of SD events. The ND distribution is normalized to the total number of SD 

events. 
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FIG. 14: The average n distribution, n* = (if ctl +rf ct2 )/2, 
of the two highest E T jets for SD and ND events. The ND 
distribution is normalized to the total number of SD events. 
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FIG. 16: The multiplicity distribution in the MPp calorimeter 
for SD and ND events. The ND distribution is normalized to 
the total number of SD events. 



data samples; in lVIBl we present ratios of SD to ND pro- 
duction rates and extract the diffractive structure func- 
tion; and in IVICi we report results on t distributions. 



A. Kinematic distributions 

The presence of a rapidity gap in diffractive events 
leads to characteristic kinematic distributions. Here, we 
compare the SD and ND distributions of mean transverse 
jet energy, mean jet pseudorapidity, jet azimuthal angle 
difference, and MiniPlug multiplicity. 

a. Transverse jet energy. The mean dijet transverse 
energy, E* T = (E 3 T etl + E 3 T et2 )/2, is presented in Fig. EE! 
for SD and ND events. The total number of ND events 
is normalized to that of the SD events and the associ- 



ated statistical uncertainties are rescaled. The SD and 
ND distributions are very similar. The slightly narrower 
width of the SD distribution is attributed to the lower 
effective collision energy in SD (JP-p collision), and the 
larger fraction of jets at 77 > in SD compared to ND 
events, which tend to have a lower average Et [27[. 

b. Jet pseudorapidity. The mean dijet pseudorapid- 
ity distribution of the two leading jets, ij* = (?7 jctl + 
?7 jct2 )/2, is shown for SD and ND events in Fig. EES The 
ND events are centered around 77* = 0, while the SD dis- 
tribution is shifted towards positive values (proton-beam 
direction) due to the boost of the center of mass system 
of the P-p collision. 

c. Azimuthal angle correlations. The distributions 
of the azimuthal angle difference between the two leading 
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jets, Acj) = |(^ ctl — (jy >ct2 \, are shown in Fig.[TS]for SD and 
ND events. The SD dijets are more back-to-back than 
the ND ones, as would be expected from the reduced 
available subenergy of the diffractive system. 

d. MiniPlug multiplicity. The MPp multiplicity dis- 
tribution for SD and ND events is shown in Fig. [16] Mul- 
tiplicities are evaluated by counting the number of peaks 
above calorimeter tower noise levels using "seed" towers 
with a minimum transverse energy of 20 MeV (3l| . Such 
peaks originate from three sources: (a) charged hadrons 
traversing the MP without interacting, (b) charged and 
neutral hadrons interacting in the material of the MP 
plates, and (c) electromagnetic energy deposited by e =t /7 
particles. The SD events have smaller multiplicities in 
comparison with the ND events, as expected from the 
reduced center of mass energy in P-p relative to p-p col- 
lisions. 



B. Ratio of single diffractive to nondiffractive 
production rates and the diffractive structure 
function 

The ratio TZ = R SD/ND (a;, Q 2 , £, t) of the SD to ND 
dijet production rates, which in leading order QCD is 
proportional to the ratio of the corresponding structure 
functions (see Sec. [TIJ) , is measured as a function of a;L, 
the Bjorken-a; of the struck parton of the antiproton, and 
Q 2 s» (E T ) 2 . For each event, x p B ^ is evaluated from the 
Et and 77 values of the jets using the formula 

3 jets 

x % = -7= s £ ( n ) 

V 2 — 1 

where the sum is carried out over the two leading jets plus 
a third jet of Et > 5 GeV, if present. Theoretically, the 
sum should be over all jets in the final state, but the frac- 
tion of events with more than three jets of Et > 5 GeV is 
relatively small and including them in the evaluation of 
x B j does not significantly affect the obtained results [I?} ■ 

Jet energies are measured using an algorithm based on 
measuring the "visible" energy deposited in the detector 
within a cone of radius R CO ne — \/ 5rf + 5<f) 2 — 0.7 and 
applying appropriate corrections [26J. Both the SD and 
ND jets contain within the jet cone an amount of un- 
derlying event (UE) energy from soft spectator partonic 
interactions, which we subtract from the measured jet 
energy. For diffractive events, where a large fraction of 
the event energy is carried away by the recoil antiproton, 
the amount of UE energy is expected to be smaller than 
in ND events. 

The UE energy in SD events was measured using the 
sample of RPS inclusive triggered data. To suppress 
overlap backgrounds, we required the events to be in 
the region of 0.01 < £^ AL < 0.14, which is dominated 
by SD events from a single pp interaction. We then se- 
lected events with only one reconstructed vertex with 



l^vtxl < 60 cm (~ 1 a) and measured the YjEt of all cen- 
tral calorimeter towers within a randomly chosen cone of 
radius i? CO nc = 0.7 in the region 0.1 < \q\ < 0.7. The UE 
energy in MB (ND dominated) events was also measured 
using the random cone technique. From these two mea- 
surements we obtained an average UE Et of 0.90 GeV 
(1.56 GeV) for SD (ND) events. 

1. XBj dependence 

The ratio TZ of the number of SD dijet events per unit 
£ over the number of ND events of the Jet5 sample is 
corrected for the effect of multiple interactions and for 
the RPS detector acceptance. 

To account for multiple interactions contributing to 
the ND sample, the ND XBj distribution is weighted by 
the factor / 1_lnt = exp[— c • C], where C is the instan- 
taneous luminosity and c = (0.34 ± 0.06) x 10~ 31 cm 2 s 
(see Sec. IV Dp . This correction is not applied to the SD 
events, since contributions from additional interactions 
shift £p into the BG region. The effect of overlaps of a 
SD and a ND event with jets in both events, which would 
tend to shift the xsj distribution towards higher values, 
is estimated to be on the 0.1% level and is neglected. 

The number of SD events is corrected for the RPS ac- 
ceptance, which is estimated from a beam optics simu- 
lation to be e = 0.80 ± 0.04 (syst) for 0.03 < £ < 0.09 
and t > 1 GeV 2 [3, and the obtained SD and ND data 
samples are normalized to their respective integrated lu- 
minosities. 
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FIG. 17: The ratio of single-diffractive (SD) to nondiffractive 
(ND) dijet-event rates as a function of XBj ■ The quoted overall 
systematic uncertainty of ±20% for the Run II ratio is due to 
the uncertainties ATZo and Ar listed in Table ITlTl The overall 
systematic uncertainty in the Run I ratio is ±25%. 

The results for the ratio TZ for events within 0.03 < £ < 
0.09 are presented in Fig. [17] as a function of log 10 XBj, 
along with the results obtained in Run I for 0.035 < £, < 
0.095 and \t\ < 1 GeV 2 [To|. The Run I results have 
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an estimated ±25% systematic uncertainty in the overall 
normalization [io| . The shapes of the two distributions 
are in good agreement, except in the high XBj region 
where the Run I distribution is seen to extend to higher 
XBj values. This behavior is expected from the difference 
in (E T ) and rj acceptance for jets in Run I (Run II), which 
is (E T ) « 7 GeV (12 GeV) and |r?| < 3.7 (|r?| < 2.7). This 
difference results in a larger x Bj reach in Run I by a factor 
of (12/7) x e 3 - 7 - 2 - 7 « 2, as estimated using Eq.QU 

A fit to all Run II data in the range 0.03 < £ < 0.09 
using the form TZ = TZq • x r B j subject to the constraint 
P = (x/£) < 0.5 Qil yields r = -0.44 ± 0.04 and TZ = 
(8.6 ± 0.8) x 10~ 3 . This result is compatible with the 
Run I result of r = -0.45 ± 0.02 and TZ = (6.1 ± 0.1) x 
10~ 3 obtained by a constrained fit of the form TZ = TZq ■ 
(xbj /0-0065) r designed to match the data at the average 
value of 2^/0.0065. This constraint accounts for the 
lower TZq value of Run I. 

The systematic uncertainties in TZq and r are 
listed in Table Mil The causes of uncertainty investi- 
gated include the underlying event, the energy scale, the 
calorimeter tower Et thresholds, overlaps in ND events, 
instantaneous luminosity C, bunch-by-bunch variations 
in C, RPS acceptance, and "splash" events. 

a. Underlying event. The underlying event energy 
is subtracted from the jet energy when the jet energy 
corrections are applied. The results presented are for a 
±30% variation of the UE energy correction, which is suf- 
ficient to cover the uncertainty for jets depositing energy 
near the interfaces or the outer edges of a calorimeter. 

b. Energy scale. The effect of the energy scale of the 
CCAL, PCAL and MP calorimeters on £? AL and thereby 
on IZo and r was estimated by changing the CCAL and 
PCAL jet energies by ±5% and the MP tower energies 
by ±30%, based on studies of inclusive jets and compar- 
isons of MP Et distributions with expectations based on 
simulations. 

c. Tower Et threshold. Tower energy threshold ef- 
fects would generally be expected to cancel out in mea- 
suring the ratio of rates. However, due to the differ- 
ent UE event contributions in SD and ND events, the 
tower thresholds applied could affect the result. The un- 
certainty in 7Z for jets with mean transverse energy of 
E T > 10 GeV (E T > 12 GeV) due to tower threshold ef- 
fects is estimated to be ±1% (±2%). The effect on both 
1Z, and on r for our total SD event sample is ±1%. 

d. Overlaps in ND events. To account for event 
overlaps in ND events occurring at high instantaneous 
luminosities, each event is weighted by a factor of w — 
exp[ -c ■ £ ], where c = (0.34 ± 0.06) x 10~ 31 cm 2 s (see 
Sec. IV D]) . The uncertainty in c has a ±8% effect on the 
determination of TZq. 

e. Instantaneous luminosity, C. We estimate an un- 
certainty of ±3% in the determination of TZo due to an 
uncertainty of ±6% in the determination of C during the 
data collection period. 

/. Bunch-by-bunch variation in C. A bunch- by- 
bunch variation of ^50% in C over the collection period 



of our data sample contributes a ±4% systematic uncer- 
tainty in TZq. 

g. RPS acceptance. Using a MC simulation, the un- 
certainty in TZq as estimated to be ±10%. 

h. Splash events. By evaluating the differences 
among data subsamples, a systematic uncertainty of ±6% 
is estimated due to the cut applied to remove the splash 
events (see Table HVf . 

Added in quadrature, the above uncertainties yield 
A1Z (syst) = ±18% and Ar (syst) = ±6%. 



TABLE III: Systematic uncertainties ATZo and Ar in TZo and 
r of the ratio of SD to ND events parametrized as TZ = TZo ■ 
x r B j. The uncertainties are estimated either by varying the 
parameter associated with the source of error by an amount 
denoted by "Variation," or in the cases marked as "N/A" (not 
applicable) by a direct calculation using data or simulated 
event distributions as explained in the text. 



Source of uncertainty 


Variation ATZo 


Ar 


Underlying event 


±30% 


±3% 


T-5% 


Central/plug cal. energy scale ±5% 


±8% 


-3% 


MiniPlug energy scale 


±30% 


T-6% 


±1% 


Tower Et threshold 


±10% 


±1% 


±1% 


Overlaps in ND events 


±20% 


±8% 


< ±1% 


Instantaneous luminosity 


±6% 


±3% 


N/A 


Bunch-by-bunch luminosity 


±50% 


±4% 


N/A 


RPS acceptance 


N/A 


±10% 


N/A 


Splash events 


N/A 


±6% 


N/A 


Total uncertainty 




±18% 


±6% 



2. Q 2 dependence 

As discussed in Sec. IIV1 special data samples collected 
with dedicated triggers are used to extend the measure- 
ment of the jet energy spectrum to E 3 T et ~100 GeV. Re- 
sults for the ratio TZ obtained from an analysis of these 
samples similar to that described above are presented in 
Fig. M for Q 2 = (E T ) 2 in the range of - 100 < Q 2 < 
10 4 GeV 2 . These distributions are mainly affected by 
the overall uncertainty in normalization, as the relative 
uncertainties among different Q 2 bins cancel out in mea- 
suring the ratio. A factor of < 2 variation among all dis- 
tributions is observed over the entire Q 2 range, as com- 
pared to a factor of ~ 10 4 over the same range between 
the individual SD and ND distributions shown in Fig. [TBI 
The results of fits performed using the form TZ = TZq ■ x r B j 
in the region 1 x 10~ 3 < xsj < 2.5 x 10 -2 , where the 
constraint j3 = (x/£) < 0.5 is satisfied, are presented 
for various jet E T intervals in Table IIVI Both TZq and r 
are constant within the uncertainties over the range of 
10 2 < Q 2 < 10 4 GeV 2 . This result indicates that the Q 2 
evolution in diffractive interactions is similar to that in 
ND interactions (42l |. 
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FIG. 18: The ratio of diffractive to nondiffractive dijet event 
rates as a function of XBj (momentum fraction of parton in the 
antiproton) for different values of Q 2 « {E T } 2 . The quoted 
overall systematic uncertainty of ±20% is due to the uncer- 
tainties A7£o and Ar listed in Table ITTT1 



TABLE IV: Fit parameters of the ratio of SD to ND produc- 
tion rates for events in different E T bins. The ratios are fitted 
to the form 72. = 7l • x Bj for 1 x 10" 3 < x B j < 2.5 x 1CT 2 . 



Jet E T [GeV] ( 


Q 2 GeV 2 


TZo 




r 


8 < E* T < 12 


100 


(8.6 ±0.8) x 10" 


-a 


-0.44 ± 0.04 


18 < Et < 25 


400 


(8.0 ± 1.6) x 10" 


-3 


-0.48 ±0.05 


35 < Et < 50 


1,600 


(6.3 ± 1.8) x 10" 


-3 


-0.60 ±0.07 


50 < E* T < 70 


3,000 


(5.5 ±5.0) x 10" 


-3 


-0.64 ±0.22 


70 < E T < 90 


6,000 


(7.0 ±7.0) x 10" 


-3 


-0.58 ±0.26 



C. t distribution 

The t distribution is measured for both soft and hard 
SD data samples using the inclusive RPStrack and the 
RPStrack-Jet5 (Jet20, Jet50) data. In Sec. IVI C 11 we 
discuss the t distribution in the region of < \t\ < 1 GeV 2 
as a function of jet (E^) 2 Q 2 and characterize the 
slope parameter at \t\ ~ as a function of Q 2 . Then, in 
Sec. IVI C~2l we examine the shape of the t distribution in 
the region up to \t\ — A GeV 2 and search for a diffraction 
minimum. Crucial for the measurement at \t\ > 4 GeV 2 
is the background subtraction. 

a. Background evaluation. The background is eval- 
uated from the data by taking advantage of the asym- 
metric vertical position of the RPS detectors during data 
taking. As shown in Fig. [SJ the dynamic alignment 
method yields Y ff sot « 0.25 cm and X g sot ~ 0. The 
nonzero value of Y ff se t results in an asymmetric t range 
in the acceptance. Since the RPS detectors are 2 cm 
wide , and given that t oc —Op, where 9p ocY ff set , the \t\- 
range of the RPS acceptance is expected to be a factor 
of [(1 + 0.25)/(l - 0.25)] 2 = 2.8 larger in the Y track > Y 



than in the Y tra ck < Yo region, where Yrj is the Y= 
horizontal centerline of the detector. 
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FIG. 19: t distributions of SD RPS tra ck dijet events within 
0.05 < £? PS < 0.08 for Y tra ck >Y (up-pointing-triangles) 
and Ytrack <Yo (down-pointing-triangles), where Yo is the 
centerline of the RPS fiber tracker. 



Figure [TO] shows t distributions of SD RPStrack data 
within 0.05 < ^ PS < 0.08. Two distributions are 
shown, one for Y trac k > Y and another for Y trac k < Y . 
The Ytrack > Yo distribution falls nearly exponentially 
from t = down to a value ~ 2.3 GeV 2 and be- 
comes consistent with an average flat background of 
iVbg ~ 20/GeV 2 in the unphysical region of \t\ > \ti\, 
in which the reconstructed track is outside the RPS ac- 
ceptance. We have verified that the event rate at \t\ > t\ 
scales with that at t ~ for runs of various instantaneous 
luminosities, and therefore conclude that these events 
represent a background associated with a pp interaction. 
However, since the reconstructed track is outside the RPS 
acceptance, its origin appears to be a particle from the 
pp interaction that suffered a secondary collision before 
reaching the RPS and produced a particle that traversed 
the RPS. Such secondary collisions would be expected to 
produce a flat distribution in the vicinity of the RPS in a 
plane perpendicular to the beam, consistent with the flat 
distribution observed at \t\ > 2.5 GeV 2 for Ytrack > Y . 

Based on the estimated value of \ti\ ~ 2.3 GeV 2 and 
the scaling factor of 2.8 for the detector acceptance due 
to the vertical misalignment, the unphysical region for 
the Yt rac k < Y data is expected to be at \t\ > \t 2 \, 
where = 2.3 x 2.8 = 6.5 GeV 2 . However, in order 
to reduce systematic uncertainties on acceptance correc- 
tions arising from (t, ^-resolution and detector edge ef- 
fects we present in Sec. lVIC~2"1 |£| distributions only up to 
\t\ = 4 GeV 2 . 

b. Corrections. The data are corrected by subtract- 
ing the background and dividing by the RPS acceptance. 
In principle, the background subtraction should be per- 
formed on an event-by-event basis, since both back- 
ground and acceptance depend on both £ and t. The 
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acceptance, A(£,t), is obtained from a beam-optics sim- 
ulation. For the data in the plateau region of the £^ PS 
distribution of Fig. [7] within 0.05 < f < 0.08, which 
are used for evaluating t distributions, the acceptance 
^4o.05<5<o.08(i) is shown in Fig. [20] For each event, the 
raw values of £ and t are obtained from the RPS tracker 
and used to fill the bins of a histogram. Since the ac- 
ceptance is fairly flat versus £ within the t region of in- 
terest, the incremental value entered into a i-bin is re- 
duced by the average background of Nb g x At events, 
where At is the bin width, and is increased by a factor of 
A7o5<£<o osW ( a zer0 vauie i s entered if the number of 
events in a bin after subtracting the background is found 
to be < 0). 

1. t distribution for \t\ < 1 GeV 2 

In this section, we discuss the event selection require- 
ments used for the measurement of the t distribution in 
the region of |t| < 1 GeV 2 , extract the slope parame- 
ter^) from fits to the data, and comment on systematic 
uncertainties. 

Event selection. To minimize the effect of migra- 
tion of events to and from adjacent £ bins caused by res- 
olution effects, events are selected in the region of \t\ < 1 
GeV 2 and 0.05 < £? PS < 0.08, where the RPS accep- 
tance is approximately flat (see Fig. [7]). To reject over- 
lap backgrounds, these events are further required to be 
within (,p AL < 0.1, where the SD events dominate (see 
Fig. [6]). The same expression that was used in the RPS 
dynamic alignment in Sec. IV Al composed of two expo- 
nential terms with slopes b\ and 62 (Eq. [7]), is used to 
fit the data and obtain the values of the slopes. Fits 
to data with only statistical uncertainties arc shown in 
Fig. [2H for various Q 2 ranges. Systematic uncertainties 
are discussed below. 

fo-slope parameters. Results for the slope parame- 
ters 61 and 62 and for the ratios b\ /b™ 1 and ^/fc™ 01 are 
presented in Table fVl The slopes 61 and 62 for all sam- 
ples are plotted in Fig. [22] No significant Q 2 dependence 
is observed from soft diffraction of < Q 2 >~ 1 GeV 2 to 
Q 2 ~ 10 4 GeV 2 . The mean values of &i and 62 over all 
data samples are 5.27 ±0.33 (GeV/cT 2 and 1.17 ±0.17 
(GeV/c) -2 , respectively. The measured slopes of the in- 
clusive sample are in agreement with theoretical expec- 
tations [431 ]. 

Systematic uncertainties. We considered the de- 
pendence of the results on the RPS fiber-tracker thresh- 
olds, the instantaneous luminosity, and the beam-store 
and run number. 

a. RPS fiber-tracker thresholds. Tracks in the RPS 
are reconstructed from hits above fiber threshold. A 
threshold value of 30 ADC counts is used for all three 
RPS stations. The slope values change by ±1% (— 1%) 
when the threshold is lowered down to 20 (raised up to 
35) ADC counts, representing a decrease (increase) that 
doubles (reduces by 50%) the number of noise hits in 



the fiber tracker. A systematic uncertainty of ±1% was 
assigned for this effect. 

b. Instantaneous luminosity. The RPS inclusive 
sample is divided into two luminosity- level subsamples: 
one consisting of events collected at low instantaneous 
luminosities, C < 1.5 • 10 31 cm _2 s _1 , and the other of 
events collected at C > 3.5 • 10 31 cm~ 2 s . A less than 
2% difference between the slopes extracted from the two 
samples is observed, which is well within the statistical 
uncertainty. 

c. Beam conditions. Variations in beam conditions 
among various beam stores and/or during various time 
intervals of a given store may affect the slope measure- 
ment. To estimate the uncertainty associated with vary- 
ing beam conditions, we measure the 6-slopes of data 
subsamples from several beam stores and from different 
periods (run number) within a store. The values for 61 
and 62 obtained were within the corresponding statisti- 
cal uncertainties. A systematic uncertainty equal to the 
statistical uncertainty was estimated. 

d. Dynamic alignment. A systematic uncertainty of 
5% is estimated from the fits shown in Fig.[5]for the slope 
61 arising from the uncertainties of ±60 /mi in the (X,Y) 
position of the RPS detectors relative to the p beam. 

The systematic uncertainties are listed in Table IVT1 



2. t distributions for \t\ < 4 GeV 2 and search for a 
diffraction minimum 

We extend the analysis to measure t distributions 
up to \t\ = 4 GeV 2 and search for a diffraction minimum 
using inclusive SD RPS tra ck and RPS tra ck-Jet5 (Jet20) 
events within 0.05 < £ PPS < 0.08. This £-range corre- 
sponds to a mean mass for the diffractively dissociated 
proton (M x ) >~ £y/s « 500 GeV. The dijets in the Jet5 
(Jet20) data sample hav e (Q 2 ) w 225 (900) GeV 2 . As 
already discussed in Sec. IVI CI the background in these 
data samples is Nt g — 20/GeV 2 , as estimated from the 
unphysical region of \t\ > \t\_\ = 2.5 GeV 2 for Y trac k > Y 
in Fig. HU 

Figure [23] shows the sum of the Y trac k > Y and 
Yt ra ck < Y distributions for the RPS inclusive and 
RPS-Jet20 data samples after background subtraction 
and acceptance corrections. The distributions are pre- 
sented in a variable-bin-width- histogram format. In in- 
crementing the histograms, a fraction is subtracted from 
each entry equal to the average background fraction in 
that bin, and the acceptance correction is then applied 
based on the RPS-measured values of £ and t. The 
dashed curve shown in the figure represents the electro- 
magnetic form factor squared of the Donnachic-Landshoff 
(DL) model [43j, Fi(t) 2 , normalized to the RPS inclusive 
data within —t < 0.5 GeV 2 . The Fi(t) form factor used 
is given by [43j]: 
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TABLE V: Slopes of t distributions for soft and hard diffractive events of the SD RPStrack data in the range 0.05 < s < 0.08 
for different (£%,) or Q 2 = (E?) 2 bins obtained from fits to the form of Eq. 0, da/dt = N ■ (Ai ■ e bl t + Ai ■ e 62 '), with 
A2/A1 = 0.11, fixed at the average value obtained in the dynamic alignment of all different event subsamples. The uncertainties 
listed are statistical. 

Event sample {E£) Q 2 bi bi bi/ 6i Ilcl 62/ & 2 ncl 

(definition) (GeV) (GeV 2 ) (GeV~ 2 ) (GeV~ 2 ) (ratio) (ratio) 
RPS incl «1 5.4 ±0.1 1.2 ±0.1 1 1 

RPSJet5 15 225 5.0 ± 0.3 1.4 ± 0.2 0.93 ± 0.08 1.12 ± 0.23 

RPSJet20 30 900 5.2 ± 0.3 1.1 ± 0.1 0.96 ± 0.07 0.93 ± 0.16 
RPS-Jet50 67 4500 5.5 ± 0.5 0.9 ± 0.2 1.00 ± 0.10 0.72 ± 0.18 
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FIG. 21: \t\p distribution shapes for SD RPStrack data of dif- 
ferent average Q 2 values within 0.05 < C^ PS < 0.08 (note the 
different £^ PS range from that in Table [I] 



TABLE VI: Systematic uncertainties in the slope parameters 
bi and 62 of the diffractive t-distributions (from Table |V)) . 



Source of uncertainty 


dbi 


Sb 2 


RPS tracker threshold 


1% 


1% 


Instantaneous luminosity 


2% 


2% 


Beam conditions 


4% 


8% 


RPS alignment 


5% 


5% 


Total (quadrature- sum) 


6.8% 


9.7% 



Fi{t) 



4m 2 , 



2.8t 



Am 2 — t 



(l-i/0.71) 
In SD, Fi(t) enters in the form [H|: 



da 



SD 



(It 



l F 1 (t) 2 exp 



2a! ■ In ~ • t 



(12) 



(13) 
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FIG. 22: The slope parameters b\ and 62 of a fit to the form 
da/dt = N ■ (A 1 -e bl ' t + A 2 -e b2 ' t ), with A 2 /A 1 = 0.11 for SD 



events of different Q 2 values (see Table [Vj ; the soft diffra 
tion (RPS inclusive) points have been placed arbitrarily at 
< Q 2 >=1 GeV 2 . 



where a' « 0.25 GeV 2 is the slope of the JP- 
trajectory 0. For \t\ < 1 (GeV/c) 2 and 0.05 < £ < 
0.08, the two-component exponential form of Eq. (J7J with 
A2IA1 = 0.11, the average value obtained in the dynamic 
alignment method, is a good approximation to that of 
Eq. (HU) with (0 = 0.065 substituted for £. 

The following features of the t distributions are no- 
table: 

• Low-t region (—t < 0.5 GeV 2 ): the RPS data are 
in good agreement with the DL curve; 

• Scale independence: the distributions of the RPS 
and RPS-Jet5 data are similar in shape; 

• High-t region (-t > 0.5 GeV 2 j: the RPS data 
lie increasingly higher than the DL curve as —t 
increases, becoming approximately flat for — t > 
2 GeV 2 . The compatibility of this observation 
with an underlying diffractive minimum at — t ~ 
2.5 GeV 2 broadened by resolution effects is dis- 
cussed below. 

The physics significance of these results is briefly dis- 
cussed below. 

a. Low-t region. The good agreement between the 
inclusive t distribution and the DL prediction in this re- 
gion serves as a basis for a search for deviations in the 
region of \t\ > 0.5 GeV 2 that could arise from a diffrac- 
tion minimum. 

b. Scale independence. The scale independence of 
the distributions supports a factorization property be- 
tween the exchange that produces the leading p and as- 
sociated rapidity gap on the one hand, and the final state 
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FIG. 23: ^distributions for two samples of SD RPStrack events 
within the region 0.05 < ^ ps < 0.08 corrected for RPS ac- 
ceptance after background subtraction: (circles) RPS inclu- 
sive and (triangles) RPS-Jet20 ((Q 2 ) ~ 900 GeV 2 ). The 
curve represents the distribution expected for soft SD in the 
DL (Donnachie-Landshoff) model j43J fEq. ll3[) normalized to 
the RPS data within -t < 0.5 GeV 2 . 



into which the proton dissociates on the other. Such be- 
havior favors models in which the hard scattering is con- 
trolled by the low- a; parton distribution function of the 
recoiling antiproton, just as in ND interactions, while a 
color-neutral soft exchange allows the antiproton to es- 
cape intact forming the rapidity gap (see, for example, 
Refs. H-S3). 

c. High-t region. In pp and pp elastic scattering at 
collider energies, a diffraction minimum (dip) in the t 
distribution is observed, with its value decreasing as y/s 
increases (see, e.g., Ref. [48j). Recently, the DO collab- 
oration reported a preliminary Tevatron Run II result 
on elastic pp scattering at y/s = 1960 GeV, in which 
a dip (broadened by resolution effects) is observed at 
—t ~ 0.7 GeV 2 followed by a maximum ("bump") at 
-t ~ 1 GeV 2 [H. A dip in the t distribution of pp or 
pp diffraction dissociation has never been reported. In 
this analysis, since the quasielastic diffractive scattering 
occurs at s' = £ • s < s, the dip, if it exists, would be 
expected to lie at a higher \t\ than in elastic scattering 
and have a Gaussian-like width due to A^-bin-size and 
t-resolution effects. 

The expected contributions to the width of a diffractive 
dip are summarized below: 



(a) £-bin width: 
ln(l/0.05) - 



from Eq. 
ln(l/0.08) 



H, A|t|c_ bil 
0.47; 



= Aln(l/0 = 



(b) 6£ and St resolutions: from Refs. [13, El, 6£ = 0.001 
and St = ±0.07 GeV 2 for (\t\) 0.05 GeV 2 with a 
dependence cx y/\t\, resulting for \t\ 2.5 GeV 2 in a 
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value larger by a factor V0.05 ~ 7 that yields 

6t Tes w 0.5 GeV 2 ; 

(c) total width expected : At|t|_2.s = A|t|g_t>i n + 5t res ~ 
1 GeV 2 , where the contributions have been added lin- 
early, since the £-bin width is not a random variable. 

These estimates show that the rather flat —t distribu- 
tions at large —t shown in Fig. [231 are compatible with a 
possible existence of an underlying diffraction minimum 
around -t ~ 2.5 GeV 2 . 



VII. CONCLUSION 

We present results from a measurement of diffractive 
dijet production in pp collisions at y/a = 1.96 TeV, p + 
p —> [jetl +jet2 + X]p + p, extracted from data collected 
in Run II of the Fermilab Tevatron pp Collider using the 
CDF II detector. We also present a measurement of the t 
distribution for single-diffractive events from an inclusive 
single diffractive data-sample up to — 1=4 GeV 2 , compare 
it to that of single diffractive dijet production, and search 
for a diffractive dip. 

A system of special forward detectors including a Ro- 
man Pot Spectrometer (RPS) for measuring the momen- 
tum of the outgoing p were used. The measured diffrac- 
tive dijet production rates as a function of x B .- (x p - 
Bjorken), Q 2 s» (Et) 2 , £p, and t p are compared with 
the corresponding nondiffractive rates at the same x p B ^ 
and Q 2 . The physics interest in these measurements is 
to elucidate the QCD nature of the diffractive exchange, 
traditionally referred to as Pomeron exchange. 

Our dijet results are extracted from data samples of 
diffractive and nondiffractive events corresponding to 310 
pb _1 of integrated luminosity. To reduce systematic un- 
certainties, the diffractive and nondiffractive data were 
collected simultaneously using the same calorimeter trig- 
ger to accept high-_Ey jets, and were similarly analyzed. 
The measurements cover the region of 0.03 < £p < 0.09, 
\t p \ < 4 GeV 2 , 0.001 < x p Bj < 0.09 and 10 2 GeV 2 < 
Q 2 < I0 4 GeV 2 . 

The E^ distribution shapes are similar for single 

diffractive and nondiffractive events, and the x B j dis- 
tribution of the ratio of single diffractive to nondiffrac- 
tive production rates is relatively flat in the region of 
/3 = x/£ < 0.5, where j3 is the fraction of the momentum 
loss carried by the parton in the Pomeron that partici- 
pates in the interaction. A fit of the form 1Z = 1Z • %Bj 
in the region of 0.001 < x Bj < 0.025, where /3 < 0.5, 
yields r = -0.44 ±0.04 and TZ Q = (8.6 ±0.8) x 10" 3 , con- 
sistent with the Run I CDF result of r = -0.45 ± 0.02 
and 1Z — (6.1 ± 0.1) x 10~ 3 within the normalization 
systematic uncertainty of ±30% common to the Run I 
and Run II data samples. 

The Q 2 dependence of the ratio of single diffractive 
to nondiffractive rates is relatively weak, varying by less 



than a factor of ~ 2 over the measured Q 2 range within 
which the individual single diffractive and nondiffractive 
distributions vary by a factor of ~ 10 4 . 

The t distributions for — t < 1 GeV 2 , for the in- 
clusive and all dijet event samples, can be fit with a 
sum of two exponential terms of the form du/dt = 
N ■ (Ai ■ e bl t + A 2 ■ e b2 '*). The slope parameters b\ and 
62 are found to be independent of Q 2 over the range of 
- 1 GeV 2 (inclusive) < Q 2 < 10 4 GeV 2 . The mean val- 
ues of 61 and 6 2 are 5.27 ± 0.33 GeV~ 2 and 1.17 ± 0.17 
GeV~ 2 , respectively. Using these slopes, the above 
double-exponential fit is in good agreement with the pre- 
diction from the Donnachie-Landshoff (DL) model [43[ 
for -t < 0.5 GeV 2 . 

The search for a diffraction minimum in the t distri- 
bution is conducted using data with 0.05 < £,p PS < 0.08 
and —t up to 4 GeV 2 . Three event-samples are studied, 
the inclusive sample and two dijet samples of (Q 2 ) ~ 225 
and w 900 GeV 2 , respectively. The distributions for all 
three samples at —t > 1.5 GeV 2 display a flattening- 
out relative to the electromagnetic form factor used in 
the DL model, possibly due to a diffraction minimum at 
—t ~ 2.5 GeV 2 filled by i-resolution and ^-dispersion ef- 
fects. At — t ~ 2 GeV 2 , the measured event-rate for the 
inclusive-event sample is a factor of ~ 10 larger than the 
DL-model prediction normalized to the data distribution 
in the region of —t < 0.5 GeV 2 . 

The relatively flat x B j distribution and the small Q 2 
dependence of the diffractive to nondiffractive ratios, 
combined with the Q 2 independence of the t distribu- 
tions, favor models of hard diffractive production in 
which the hard scattering is controlled by the parton- 
distribution-function of the recoil antiproton while the 
rapidity-gap formation is governed by a color-neutral soft 
exchange [44l]-[47l|. 
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